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Abstract: Marginal likelihoods for the cosmic expansion rates are evaluated using the 
recent 'Constitution' data of 397 supernovas, thereby updating the results in some previous 
works. Even when beginning with a very strong prior probability that favors an accelerated 
expansion, we end up with a marginal likelihood for the deceleration parameter qo peaked 
around zero in the spatially flat case. This is in agreement with some other analysis of the 
Constitution data. It is also found that the new data significantly constrains the cosmic 
expansion rates, when compared to the previous analyses. Here again we adopt the model- 
independent approach in which the scale factor is expanded into a Taylor series in time 
about the present epoch; for practical purposes, it is truncated to polynomials of various 
orders, in different trials. Though one cannot regard the polynomials thus obtained as 
models, in this paper we evaluate the total likelihoods (Bayesian evidences) for them to find 
the order of the polynomial having the largest likelihood. Analysis using the Constitution 
data shows that the largest likelihood occurs for the fourth order polynomial and is of 
value ps 0.77 x 10" 102 . It is argued that this value, which we call the likelihood for the 
model-independent approach, may be used to calibrate the performance of realistic models. 
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1. Introduction 

A decade ago, the apparent magnitude-redshift (m — z) data of type la supernova (SN la) 
were declared to be indicating that our universe is expanding with an acceleration |||. 
This discovery is considered to be the most important and startling one in cosmology, after 
the 1929 observation by Hubble that the universe is expanding, and the 1964 discovery 
by Penzias and Wilson that there is a cosmic microwave background. The discovery of 
accelerated expansion is startling for it demands the presence of substantial amount of 
some unknown 'dark energy' in the universe, with repulsive pressure. For the analysis of 
the data, both the supernova search teams which first reported this observation assumed the 
validity of the LCDM model in cosmology, in which the universe contains a cosmological 
constant A or an equivalent vacuum energy p\ (with an equation of state p\ = —pa), 
along with nonrelativistic matter. The latter includes some unknown cold 'dark matter' 
too. Analyses were also performed under the speculation that the dark energy has the 
equation of state pd e = wpd e , with it; as a time- varying quantity. The Chevallier-Polarski- 
Linder (CPL) ansatz j|, |I| is the popular one in this connection. Later it was speculated 
that there is interaction between the energy components |(|. It may even be asked 
whether there is only one component in the 'dark energy' at the present epoch 

We may note that in all the above analyses, the following assumptions are made: 
(1) The Einstein equation in general relativity is valid for the cosmos (which are the 
Friedmann equations). (2) The universe is homogeneous and isotropic at very large scales 
and hence spacetime possesses a Robertson- Walker (RW) metric. (3) There can be one 
or more components in the cosmic fluid and in the latter case, the components may or 
may not be separately conserved. (However, the total energy density will be conserved.) 
Specific assumptions in this regard, on each component and their equations of state, are 
prerequisites to begin the analysis, for we need a solution a(t) of the Einstein equation in 
this case (where a(t) is the scale factor of expansion) at our disposal. Such assumptions 
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and the resulting solution a(t) constitute a 'cosmological model'. It shall be noted that a(t) 
often contains some free parameters, but the specific values assumed by these parameters 
alone do not lead to different models, unless those values are predictions in them, on the 
basis of some fundamental principles. 

When we make such model-based analyses, one of the tasks is to find the model which 
enjoys the largest support from the data. The merit of a new model is adjudged by compar- 
ing it with any existing ones, using the data related to the phenomena under consideration. 
For this, the conventional approach is to perform a maximum likelihood ratio test where 
one has to find the lowest \ 2 f° r the models concerned. But it is now recognized that 
the use of Bayes theory is the more reasonable approach in problems such as cosmology 
g |, 0]. In Bayesian model comparison, we take two models as rival hypotheses (say, 
Mi and Mj) and try to compare them by evaluating odds ratios between their posterior 
(i.e., after analyzing the data) probabilities, given the data D and also some background 
information /. To obtain these ratios, the Bayes's theorem is made use of. In the cosmo- 
logical literature, a large amount of work which use Bayes theorem is now reported [See 
12, 13] for some reviews]. But a clear distinction between Bayesian parameter estimation 



and Bayesian model comparison is not found in most of the reported works. It shall here 
be noted that cosmological models which are truly distinct (in the sense of the definition 
of the term given above) were first compared using Bayesian theory in [ 10 1 . 

We note that since there can be infinitely many solutions a(t) based on our specific 
assumptions regarding energy densities, it is not possible to evaluate any one model's pos- 
terior probability. Instead, as mentioned above, one can perform a model comparison. 
Comparison is possible because here we need only to find the ratios between the probabil- 
ities for obtaining the data D in the various models, multiplied by any prior odds. The 
probability for obtaining the data D in a model is often referred to as the 'likelihood for 
the model' or the 'Bayesian evidence'. 

On another front, it was pointed out that the most appropriate way to measure the 
acceleration of the universe is to resort to a model-independent approach. In the con- 
ventional model-based analyses of m — z data of SN la, the accelerated expansion of the 
universe was an indirect inference based on the best fit values of parameters, such as the 
density parameters f2 m and J1a in the LCDM model. In the model-independent case, the 
scale factor a(t) is expanded as a Taylor series in time about the present epoch |l4], and 
the marginal likelihoods of its coefficients are computed using the data. But practically, 
we have to truncate the series to some finite order and hence our basic assumption is that 
a(t) is expressible as a truncated Taylor series or polynomial. Evaluating the deceleration 
parameter by adopting this method, it was confirmed model-independently using the SN 
data that the universe is undergoing an accelerated expansion [|14|]. Clearly, this Taylor 
expanded scale factor is not a realistic cosmological model, since here homogeneity and 
isotropy are the only assumptions made for the universe. 

Model-independent approaches which use Taylor expansion in terms of redshift z have 
also gained wide attention in recent years [See for eg. |l|, H, A great advantage 
of the present approach of expanding the scale factor in terms of t about the present epoch 
to, when compared to expansion in z about z = is that the former converges for all times, 
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whereas the latter converges only for | z |< 1 [17]. In the present paper, even the lookback 
time T = t — to is evaluated by numerically solving Eq. (13) in Ref. [15], which involves 
a Taylor series in time; here we do not use an expansion in z at all. Hence there is no 
convergence problem [17, 18 1 in the present work. 

It was mentioned above that in Bayesian model comparison one evaluates only the 
relative merits of models on the basis of data and it is not possible to compute the posterior 
probability for a model. Also, Bayesian model comparison does not give the best possible 
model, for it only compares the available ones. In this circumstance we ask whether one 
can find some standard, which can be considered as a minimum requirement, for a model 
to be termed successful. In the lowest x 2 P er degree of freedom approach, such a crude 
standard exists; a model is considered a reasonable good fit if the x 2 per degree of freedom 
is less than or nearly equal to unity. Several analyses of SN data in cosmology still use 
this standard [for eg. [20, 21 1]. We attempt in this paper to set a similar standard in 
Bayesian analysis. Obviously, the relevant quantity in this context is the likelihood for 
the model or the Bayesian evidence. We propose to evaluate this quantity for the scale 
factor a(t) in its Taylor series form about the present epoch, thus combining the Bayesian 
and model-independent approaches. Eventhough one cannot regard the truncated Taylor 
series form of the scale factor as a model, the likelihood for polynomials of various orders 
may be obtained, and thereby one can obtain the order of the polynomial with the largest 
likelihood. That such a maximum exists can be seen from the fact that the Bayesian 
approach has a built-in mechanism to implement the Occam's razor; i.e., it favors simpler 
models when compared to complicated ones with more parameters, unless the latter shows 
significantly better performance. We have assumed that the scale factor of the universe has 
a unique Taylor expansion with definite values for its coefficients, which we attempt to find 
with the help of data. Varying the coefficients in the series arbitrarily will certainly affect 
X 2 ; they are not unconstrained parameters. When the order of the truncated Taylor series 
becomes very large, the number of coefficients in it too becomes large so that Occam's razor 
forces the likelihood of the polynomial to tend to zero. Therefore it should be possible to 
find the order of the polynomial that maximizes the likelihood. We here argue that this 
maximum value of the likelihood/evidence for the Taylor expansion, which may be termed 
the likelihood for the model-independent approach, can be set as a standard for model 
comparison. 

It shall here be noted that a great practical use of scientific theories is that the equations 
they provide save us from keeping large amount of raw data, for use in future applications. 
Therefore it is natural to expect that the likelihood for a fundamental scientific theory 
exceed that of a truncated Taylor series of the unknown function or at least be equal to 
that of such a series. In other words, only those realistic models, which have likelihoods 
greater than or at least equal to that in the model-independent approach can be considered 
as successful. Even in cases where it is not possible to find the order of the polynomial 
that maximises the likelihood for some practical reason, it is reasonable to demand that 
realistic models perform better than each of the (low order) polynomials we have worked 
out. 

Bayesian model-independent approaches of the kind proposed here are pursued by 
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some authors, though not in a systematic way as in [15] or as proposed to do in this paper. 
For instance, Guimaraes, Cunha and Lima [18] have compared the realistic flat LCDM 
model with simple kinematic models, such as those based on three simple parametrisa- 
tions for the deceleration parameter etc., using the 307 SN la Union compilation data 
set. They found that even very simple kinematic models are equally good to describe this 
data, when compared to LCDM model. We develop this procedure further and attempt to 
achieve a systematic calibration of cosmological models, by first evaluating the likelihoods 
for polynomials of various orders and then finding the largest possible value of the likeli- 
hood for a Taylor expansion. For calibration purpose, we compare other model likelihoods 
with this value using the Bayesian method. Moreover, in this work we use a more recent 
'Constitution' supernova data |2j| which contains 397 objects. 

Another important work we report in this paper is that of updating the marginal likeli- 
hood for each of the expansion coefficients found in [14, 15]. This is performed for the case 
of a fifth order polynomial. The new marginal likelihoods for its coefficients give valuable 
information regarding the expansion history of the universe. An interesting result in this 
connection is that even when beginning with a very strong prior probability that favors an 
accelerated expansion, we end up with a marginal likelihood for the deceleration parameter 
qo peaked around zero in the spatially flat case. This result is in agreement with some other 
analyses [See for eg., [20|]. It is also found that the new data significantly constrains the 
cosmic expansion rates appearing in the Taylor expansion, when compared to the previous 
analyses. We also note that successive terms in the series decrease sufficiently fast, thereby 
verifying the assumption of Taylor expansion. It is expected that in the near future, as 
the SN dataset becomes large enough, these coefficients get sharply peaked marginal likeli- 
hoods and become the most basic model-independent description of the expansion history 
of the universe. 



2. Bayesian model-independent approach 

For parameterized models with parameters a, (3, .., the likelihood for the model Mj or the 
Bayesian evidence [denoted as £(Mj)], which is the probability for the data D given the 
truth of the model Mi, can be evaluated as 



p{D\M u I) = C{Mi) 
= J da Jdj3...p(a,P,...\M i )£ i (a, /?,...), (2.1) 

where p(a, (3, ...\Mi) is the prior probability for the set of parameter values a, (3, .. given the 
truth of Mi and £i(a, (3, ■■■) is their likelihood function. The likelihood function is often 
taken to be Q 

d(a, 0, ...) = exp [- X f (a, (3, ..)/2] . (2.2) 

where xf is the X 2 ~ s tatistic. Another quantity of interest is the marginal likelihood for 
any one parameter, say a, in a model. This is obtained by integrating the integrand in 
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( |2.1| ) over all parameters, except a. Thus in the above case, the marginal likelihood for the 
parameter a can be obtained as 

A(a) = J dp J d 1 ...p{a,P, 1 ...\M i )C i {a,P, 1 ...). (2-3) 

While using the apparent magnitude-redshift data, the observable is the distance mod- 
ulus n = 5 log (D/lMpc) + 25, a function of redshift z and contains parameters a, /3, etc. in 
the model concerned. D(z; a, /?..) is called the luminosity distance and D/l Mpc indicates 
that it is expressed in units of megaparsec. 

In the model-independent analysis of SN data |l4|, the scale factor of the universe 
is expanded into a Taylor series in t about the present epoch to- With t — to = T, where to 
is the present time, the Taylor series can be written as 



a(t + T) = a x (2.4) 
1 + HqT - ?°Jp-T 2 + ^T 3 - ^T 4 + ^pT 5 + ..' 

Note that T assumes negative values. Here the parameters in the theory are the present 
value of the scale factor ao, the Hubble parameter Ho = lOOh km s _1 Mpc— 1 , the deceler- 
ation parameter q^, higher order expansion rates such as tq, sq, uq, etc. and the curvature 
scalar k = dbl. (The spatially flat k = case is incorporated by including sufficiently 
large values of ao-) One of our tasks is to deduce the values of these parameters from the 
observational data. 

For a light pulse emitted from an SN situated at the coordinate r% at time t\ and 
reaching us at r = at time to, the RW metric allows one to write 



k o(0 Jr, (l-kryrr 
For a k = RW metric, this can be used to obtain 

= cdt_ f° dT 
ri J tl a(t) J Ti a(to + TY { ■ ] 

Similar expressions can be found for k = ±1 cases too. With this, we may compute the 
luminosity distance D = riao(l + z). An important part of the calculation is the solution 
of the following equation, used to find T\ in terms of z, for each combination of parameter 
values. This is done in a direct and purely numerical way: 

1 + ^ = ^%T (2-7) 

We may thus obtain the distance modulus \x = 5 log (D/lMpc) + 25. Here D and hence \x 
are functions of z and contain parameters k, ao, H , q , r , so, u , etc. 

The likelihood function is now C = exp[— x 2 (k, h, ao, qo, fo, so, iio)/2] where x 2 is given 

by 
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^2 = Sfc ^ Afc - Vk(zk;k,h,ao,q ,r ..) ^ ^ ^ 

Here /}& is the measured value of the distance modulus of the k th supernova, [ik{ z k', k, h, ao,qo, ••) 
is its expected value (from theory) and is the uncertainty in the measurement. 

The likelihood for the truncated Taylor series form of scale factor can be found using 
equation (|2.1|) as 



p(h)p(a )p(qo)p(ro)p(s )p{u ) e - * 2 / 2 . (2.9) 

where p(h)p(ao)p(qo)p(ro)p(so)p(uo) is a product of Gaussian probability distributions eval- 
uated using the mean values and standard deviation of the marginal likelihoods obtained 



in previous analyses [[15] in our case] for each of the parameters. This is an approximation 

to p(h, ao, qo,ro | Mj), the prior probability to be introduced in equation (|2~l|). 

In the above, we have kept terms up to fifth order in the Taylor expansion. Increasing 
the number of terms by unity will enhance the computation time by more than an order 
of magnitude. We have performed this computation with various orders in the truncated 
series, starting with second order in which only the first three terms are kept. For orders 
different from 5, necessary changes are to be made in the above expressions. 

3. Marginal likelihoods for the cosmic expansion rates 

We first obtain the marginal likelihoods for the various expansion rates in the case of 
truncated Taylor series of order 5. This is a repetition of the calculation in 0, [jq| , using 
the Constitution data. The marginal likelihoods we obtain here give valuable information 
regarding the expansion history of the universe. 

An important step made in the present computation of marginal likelihood is that while 



using (2.9), the marginal likelihoods in the previous analysis JlJ, 15] are taken as the prior 
probability distributions, for the corresponding coefficients. References [14, 15] have used 
flat priors, since there were no other previous work evaluating these marginal likelihoods. 
But there itself it was proposed that the posterior marginal likelihoods obtained shall 
be used for subsequent analysis and the present work is the appropriate place to make 
use of this. However, it would not be computationally feasible to use the posterior in 
the previous analysis as prior in terms of a table of values; instead, as stated above, we 
approximate those distributions by Gaussian functions with the corresponding mean and 
standard deviations obtained in |i~5f ]. A comparison with the actual plots show that this 
is a reasonable approximation for most coefficients and at any rate is a better option than 
flat priors. The product of such individual priors is the combined prior, used in equation 

(ED. 

In the present model-independent analysis, the 'Constitution' data [^] of 397 SN were 
used. We have computed the marginal likelihoods of four important expansion rates of 
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the present universe, namely qo, ro, so and uq, and the results are shown in Figs. (l)-(4). 
We have kept terms up to fifth order in this computation, but only the flat (k = 0) case 
is considered here. This is equivalent to assuming a (^-function prior for the flat spatial 
geometry. The joint prior probability used for other parameters was, as described above, 
the product of individual Gaussian functions in each parameter with mean and standard 
deviations as follows: h = 0.68 ± 0.06, q = -0.90 ± 0.65, r = 2.7 ± 6.7, s = 36.5 ± 52.9, 
and uq = 142. 7 ±320 pf. In each case, the integrations were performed in the 2a range of 
each of the parameters. We have performed variation with respect to h, though marginal 
likelihood for this parameter was not drawn. The step sizes chosen for these parameters 
were Ah = 0.01, Aq = 0.1, Ar = 1, As = 20 and Au = 100. 




Figure 1: Marginal likelihood for the parameter go (in units of 10 105 ), while using the polynomial 
of order 5. 

The results show that there is significant constraining of the parameters while using the 



new and refined data, compared to the corresponding results in [15]. It is to be reminded 
that the marginal likelihoods are not precisely probability distributions for the parameters; 
instead, they are the probability for the data, given the model and the parameter values. 
However, we here compute mean and standard deviations considering them as distributions. 
The new mean and standard deviations are the following: go = 0.04±0.30, ro = — 4.5±4.6, 
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Figure 2: Marginal likelihood for the parameter tq (in units of 10 105 ), while using the polynomial 
of order 5. 



sq = —42.8 ± 52.5, and uq = 320.5 =L 213.0. The marginal likelihood for qo obtained in 



flq l [which was not much different from that of [Q, both using the same data in [23|], is 
reproduced here in Fig. (|5|) for comparison with the distribution in Fig. (|l|). It can be seen 
that the standard deviations of each of these parameters except that of so have decreased 
substantially, which leads to our above assertion. 

It shall be noted that even when beginning with a prior centred around qo = —0.9, 
which is strongly in favor of an accelerated expansion, we ended up with this marginal 
likelihood peaked around qo « slightly towards the positive side. This casts doubt on 



the paradigm of cosmic acceleration itself. Where as the data in [23| validated the claim 
of accelerated expansion |ll|, |15]| , the new extended SN dataset in [p^ j indicates that the 
universe is neither accelerating nor decelerating. This agrees with the analysis in [^] , which 
finds that a 'coasting' (qo = 0) evolution for the universe is equally plausible. However, 
the presence of substantial amount of dark energy and dark matter would still be required 
to explain the data. 

The considerable spread in the marginal likelihoods shows that there is still enough 
freedom in choosing the values of those parameters for a best fit. In other words, even now 
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Figure 3: Marginal likelihood for the parameter so (in units of 10 105 ), while using the polynomial 
of order 5. 

there is a sizable volume in the parameter space that can have the same low \ 2 - But this 
should not be viewed as a drawback of the analysis; instead, this simply reflects the fact 
that the data are not yet accurate enough. Several recent analyses of Constitution SN data 



have reported that such freedom exists [20, 21]. This freedom in SN data was noted earlier 



in fl"4| , O], which highlights the strength of the Bayesian model-independent approach. 

Based on the mean values obtained for these parameters, we compute the successive 
terms in the series (|2.5|). With time in units of 10 17 s, the series can be written as 



1 + 2.106 x lO^T - 2.22 x 10~ 2 g T 2 + 1.55 x 10~ 3 r T 3 

-0.819 x 10~ 4 s T 4 + 3.45 x 10~ 6 u T 5 + (3.1) 

where we have taken h = 0.65. With the values of the parameter in the ranges obtained in 
the analysis, this series appears to converge even for | T | as large as ~ 3 x 10 17 s. However, 
this feature is not essential for our analysis, for we have assumed only a polynomial form 
for the scale factor. The situation was not different in [[i~J|, |l5| either. 



II II II I I II 




I I I I I I I I I I I 

f)00 o r>oo ■ noo ir>oo 



Figure 4: Marginal likelihood for the parameter uq (in units of 10 107 ), while using the polynomial 
of order 5. 

4. Likelihoods for the models and the model-independent approach 

Now we evaluate the likelihoods for various models and also for the model-independent 
approach. The realistic models considered are the flat LCDM model with constant equation 
of state parameter w = —1 and the more general flat dark energy model with the CPL 
equation of state. For the model-independent case, we have used truncated serieses of 
various orders, starting with order 2 to the largest order 6 for the scale factor. We have not 
restricted ourselves to the flat case in this analysis (whereas in the above section, all the 
marginal likelihoods were computed by assuming a fifth order polynomial and with k = 0). 

These likelihoods in different models are important since the Bayes factor is the ratio 
between them. We consider that there is positive evidence for a model, when compared 
with another one, only if its likelihood is greater than 3 times the value of the latter B. 
The evidence is considered to be strong only if its likelihood is larger than 20 times and 
very strong if it is larger than 150 times Q that of the other. Therefore it is highly desirable 
to find the likelihoods or Bayesian evidences for models and we now compute them for the 
cases mentioned above. 
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Polynomial of order 2 (up to qo) 


350 x 10 


-105 


Polynomial of order 3 (up to ro) 


90 x 10 _ 


105 


Polynomial of order 4 (up to sq) 


770 x 10 


-105 


Polynomial of order 5 (up to uq) 


30 x 10 _ 


105 


Polynomial of order 6 (up to vq) 


10 x 10" 


105 


Flat LCDM model 


9.8 x 10" 


-105 


Flat dark energy model 


0.89 x 10" 105 


with CPL equation of state 







Table 1: Likclihoods/Bayesian evidence for various mod- 
els. We have taken k = 0, ±1 in the case of polynomial 
approximations but k = for the two realistic models,. 



The likelihood for polynomials of various orders we have found are given in Table 1. It 
is seen that in this model-independent approach, the polynomial of order 4 has the largest 
likelihood ~ 770 x 10- 105 and as discussed in the introduction, this value shall be taken as 
the likelihood for the model- independent approach. 

For flat LCDM model, the like- 
lihood was 9.8 x 10" 105 and for flat 
dark energy model with CPL equa- 
tion of state, it was 0.89 x 10- 105 . 
These values, when compared to the 
likelihood 770 x 10- 105 of the model- 
independent approach, show that 
the performaces of those realistic 
models are not as good as the poly- 
nomials. 

The priors used for the coef- 
ficients in the polynomial form of 
the scale factor are the same as the 
Gaussian priors used in the previous section. For vq, we took a flat prior in the range 
—3000 < Vq < 3000 in the same manner as we chose flat priors for other parameters 
in the previous work; i.e., flat priors only for the contributing regions of the parameter 
concerned. For the realistic models, we took k = and the prior for h was the same as 
that in the previous section. Q m had Gaussian prior, with mean and standard deviations 
f2 = 0.29 ± 0.025. Integrations were performed over 2a range of these parameter values. 

The two additional parameters in the model with CPL equation of state were chosen 
to have flat priors with —1.2 < wq < 0, —12 < w\ < 0, which are judged from some recent 
computations |21|| . Variations in these priors can change the likelihoods for this model to 
some extent, which can affect our conclusions regarding the dark energy model in a similar 
manner. But we note that the lowest value for x 2 claimed by [21| for this model is 461.254 
at Q m = 0.453, wq = —0.207 and w\ = —11.316. Even if we set <5-function priors at these 
parameter values, the likelihood can only reach the maximum possible value = 6.9 x 10" 101 , 
which corresponds to this x 2 value. But the Gaussian prior peaked at Cl m = 0.29 we must 
use is far from such a 5- function distribution located at £l m = 0.453, so that there is no 
chance for this model getting a high value for the likelihood. Moreover, using the same 
(5-function distributions as priors for these parameters, while analyzing some future data, 
may turn out to be extremely harmful for this model. Therefore we resort to using the flat 
priors for wq and w±, as mentioned above, but at the same time caution that the present 
conclusion regarding the model with CPL equation of state is liable to change since we do 
not have fiducial priors for these parameters. Only after several repetitions can we reach 
such priors for them. However, the low value of likelihood obtained for the LCDM model 
is more or less robust. 

The step sizes chosen for the parameters were Ah = 0.01, Ago = 0.1, Atq = 1, Aso = 
20 and Auq = 100, as in the previous case. For other parameters, we took Avq = 1000 
and Aao = 1000. The flat case is incorporated here by including large values of qq, say up 
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to ao = 8000 Mpc and considering only k = ±1. For wo and w\, we chose Awq = 0.1 and 
Au>i = 1. 



5. Conclusion 



We assumed that a Taylor series form for the unknown variable a(t), which describes the 
cosmic expansion history, is valid and first attempted to find the coefficients in this expan- 
sion using the recent Constitution SN data. It is found that there is significant constraining 
of these parameters when compared to a previous analyses using the data in [23]. The new 
marginal likelihoods for various coefficients evaluated with the Constitution data lead us 
to expect that when more refined and abundant SN data set becomes available in the near 
future, the curves will get sharply peaked and this method of direct determination of the 
cosmic expansion history shall prove to be indispensable. 

One of the notable results obtained from the present analysis is the shift in the com- 
puted mean value of the deceleration parameter qo, from that found in the previous analysis. 
Even when we start with a prior probability distribution that strongly favors an acceler- 
ating universe, the present analysis using Constitution data provide a marginal likelihood 
peaked around the zero of the deceleration parameter. This result is in agreement with 
the analysis in (^0|. However, we reiterate that the considerable spread still found in the 
likelihoods of these parameters indicate freedom in the choice of their numerical values. 
This also is in agreement with several other analyses of the same data using alternative 
methods ||, 

Another attempt we made in this paper is the comparison of the performance of some 
realistic cosmological models with that of the model-independent approach, in explaining 
the SN data. We note that a great practical use of scientific theories is that the equations 
they provide save us from keeping large amount of raw data, for use in future applications. 
Therefore it is natural to expect that the likelihood for a fundamental scientific theory 
exceed that of a truncated Taylor series of the unknown function (the model-independent 
approach, as we refer to in this paper) or at least be equal to that of such a series. We 
here find that the two popular realistic cosmological models analyzed are much behind 
a Taylor expansion for the scale factor in explaining the data, since their Bayesian evi- 
dences/likelihoods are less than that in the latter case. We have used fiducial Gaussian 
priors for Q m in the flat LCDM model, but caution that the likelihood for the dark energy 
model with CPL equation of state is liable to change since such priors are not available for 
the parameters wq and w\. 

A distinguishing feature of our analysis is that the marginal likelihoods for each pa- 
rameter obtained in the previous case is chosen as the prior probability distribution in the 
present one, thereby implementing the Bayesian method in true spirits. The work is also 
intended as a demonstration of this fundamental requirement in Bayesian analysis. We 
have also verified that successive terms in the series decreases fast enough, justifying the 
Taylor expansion hypothesis, though this is not essential since our basic assumption was 
that the scale factor is expressible as a polynomial. 
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What we envisage here is a combined Bayesian model-independent approach. The 
application of this method answers one of the criticisms raised against Bayesian model 
comparison - that it is soft towards models with poor explanatory power, since it only 
compares the available ones. We now have a new quantity, which we call likelihood for 
the model- independent approach. Evaluation of this likelihood/evidence is only a logical 
extension of the evaluation of marginal likelihoods of those coefficients in the expansion 
and it is now demonstrated that this quantity can be used to calibrate the performace of 
cosmological models. 
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